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INTRODUCTION 



Recently, several high-precision experiments have been performed on the bound-electron 
5f-factor in H-like carbon and oxygen by the Mainz-GSI collaboration The value 

actually measured in the experiment is {gm/M), where m is the electron mass, M is the ion 
mass, and g is the (7-factor of the electron. The relative accuracy of the best experimental 
determination of this value is 5 x 10~^°, which is 4 times better than that of the accepted 
value for the electron mass |3|. Further progress is anticipated from the experimental side, 
as well as extension of measurements to the higher-Z region 0|. 

The spectacular experimental results have triggeredereat interest to the theoretical de- 
scription of the (yf-factor of a bound electron 

lisl . Ilfll 1^ . Combining experimental values with accurate theoretical predictions for 
the bound-electron gf-factor resulted in an independent determination of the electron mass 
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The current accuracy of this determination is 4 times better than that of the 



accepted value for the electron mass [3|. For the latest compilation of various contributions 







2l| for 



to the bound-electron (^-factor we refer the reader to |2Q] for H-like ions and to 
Li-like ions. 

In the present work we give a detailed description of our calculation of the self-energy 
correction to the bound-electron (yf-factor for Is and 2s states of a H-like ion. The first results 
of this calculation for the Is state were previously published in [ij], where they were used 
for the determination of the electron mass. In this paper, we extend our consideration to a 
higher-Z region and perform calculations also for the 2s state, having in mind the planned 
extension of the experiments to Li-like systems. 

Our calculation is carried out in the Feynman gauge. The relativistic units {h = c = 1) 
and the Heaviside charge units (a = e'^/in, e < 0) are used throughout the paper. We also 
use the notations = Pfj,'j^ and p = p/|p|. 

II. BASIC FORMULAS 

In this paper we will consider a bound electron in an s state of a H-like ions with a 
spinless nucleus interacting with a static homogeneous magnetic field. The bound-electron 
(yf-factor is defined by 

where /x is the operator of the magnetic moment of the electron, fiQ = |e|/(2m) is the Bohr 
magneton, ja is the total angular momentum of the electron, and is its projection. The 
lowest-order value for the q factor can be found by a simple relativistic calculation based on 

n 

the Dirac equation [22]. For an s state and the point nucleus, it yields 

+ (2) 



3 \ m ^ 
where Ea is the energy of the electron state. 

Various contributions to the g-iactor are related to the corresponding corrections to the 
energy shift, as given by 

AE = AgfXoBma, (3) 

where the magnetic field B is assumed to be directed along the z axis. In the present in- 
vestigation we are interested in the self-energy correction to the bound-electron gf-factor. It 



is diagrammatically depicted in Fig. ^ Formal expressions for the corresponding contribu- 



tions can easily be derived, e.g., by the two-time Green function method 



23|. They can be 



also obtained by considering the first-ord er p erturbation of the self-energy correction by a 



perturbing magnetic potential (5V^(x) 
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The contribution of the diagram (b) is referred to as the vertex ("ver") term and is given 

by 
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where 5V{^ = — ea- A(x), A is the vector potential A(x) = (1/2)[B xx], m = 1 — iO where 
the small imaginary addition preserves the correct circumvention of poles of the electron 
propagators, I{uj) is the operator of the electron-electron interaction, 

I{lu) = e^a^a^'D^^iu) , (5) 

= (1, ex.) are the Dirac matrices, and Df^^, is the photon propagator. The contribution 
of the diagrams (a) and (c) is conveniently divided into 2 parts that are referred to as 
the irreducible and the reducible one. The reducible ("red") contribution is defined as a 
part in which the intermediate states in the spectral decomposition of the middle electron 
propagator (between the self-energy loop and the magnetic interaction) coincide with the 
initial valence state. The irreducible ("ir") part is given by the remainder. It can be written 
in terms of non-diagonal matrix elements of the self-energy function as 

AEi, = (5a|7°S(£,)|a) + (a|7°S(£,)|5a) , (6) 

where the perturbed wave function is given by 

|5a)='f5''Mi!^^, (7) 

n 

S(£:) = S(£:) — 5m, ^(e) is the self-energy function defined as 

/oo 
duj 
'OO 

xG(£: - u;,xi,X2)a,,-D^''(u;,Xi2) , (8) 



6m is the mass counterterm, and G is the Dirac-Coulomb Green function, G{e) = (e — H) 
where H is the Dirac-Coulomb Hamiltonian. The expression for the reducible part reads 



-1 



\a) . (9) 



III. DETAILED ANALYSIS 

The irreducible contribution (jH)) is written in terms of non-diagonal matrix elements of 
the self-energy function. Its renormalization is well known and we do not discuss it here. 
For the point nuclear model, the perturbed wave function \6a) can be evaluated analytically 
in a closed form by using the method of generalized virial relations for the Dirac equation 
The explicit expression for the \6a) function can be found in 

Both the vertex and reducible contributions are ultraviolet (UV) divergent. In order 
to covariantly renormalize them, we expand the bound-electron propagators in Eqs. (jU 
and in terms of the interaction with the binding field. According to the number of 
interactions (with the nuclear Coulomb field), the corresponding contributions are referred 
to as the zero-potential, one-potential, and many-potential terms. For the vertex diagram, 
this decomposition is schematically presented in Fig. |21 We thus represent the vertex and 
reducible parts as 

AE^er = A42 + AEW + AEilp , (10) 
AE,ed = AElll + AEill + AE^'J^ , (11) 

where the superscript corresponds to the number of interactions with the binding field. It 
is often convenient to consider the vertex and the reducible terms together. In this case, we 
will use the following notation 

AE« = AE^l + AeS, , (12) 

where i = {0, 1,2+}. 

By elementary power-counting arguments one can show that all terms containing one 
or more interactions with the binding field are UV finite. Despite the fact that the one- 

n 

potential term is finite, we prefer to consider it separately, as was first proposed in p]. The 
reason for such treatment is that, when evaluated in coordinate space, the one-potential term 
yields a slowly-converging partial-wave expansion, which represents the main computational 



difficulty in the \ow-Z region. The numerical scheme developed for the computation of this 
term in j^Q] allowed to extend the summation well beyond 100 partial waves, which explains 
much better accuracy obtained in that work as compared to the ffist evaluation [5]. In our 
approach, we evaluate the one-potential term directly in momentum space without utilizing 
the partial-wave expansion. In this way we eliminate the uncertainty due to the estimation 
of the uncalculated tail of the series. 



A. Zero-potential contribution 

The expression for the zero-potential vertex term can be obtained from Eq. by replac- 
ing all bound-electron propagators by the free propagators. Writing it in momentum space, 
we have 

A^(S = -e I ^^,(p)r^(p,p')-A(p-p')V^a(p'), (13) 

where ip = V'W^j P ^i^nd p' are 4- vectors with a fixed time component, p = (e^, p), p' = {Sa, p'), 
A(p — p') is the vector potential A(x) in momentum space, 

A(p - pO = -'-i2nf [B X V^,5%p - p')] , (14) 

and r^(p, p') is the UV-finite part of the free vertex function introduced in Appendix 1X1 
We note that right from the beginning we are working with the renormalized expressions; 
the cancellation of UV-divergent terms in the sum of the zero-potential vertex and reducible 
terms can be checked explicitly. 

Substituting Eq. (fT^ into Eq. (fT!^ and separating the contribution to the g-iactor, we 
obtain 

A.iS = -2.m I |^?.(P) 

X [Vp,53(p-p') X r«(p,p')],^a(p'), (15) 

where the gradient Vp' acts on the 6 function only and the angular-momentum projection 
of the initial state is assumed to be rria = 1/2. The expression above is transformed by inte- 
grating by parts and performing the integration that involves the 6 function. The integration 
by parts divides the whole contribution into two pieces 

A^?(°) = Agi%, + Agil, (16) 



that correspond to the gradient acting on the vertex and the wave function, respectively, 

A^S,i = 2^"^ / ^a(p) m Up) , (17) 

Agi%, = -2im j ^^a(p) [Tr{p,p) x Vp]^^,(p) , (18) 

where 

s(p)= [Vp, xr(p,p')LU,- (19) 

We start with the first contribution. The function S can be expressed as 



Transformation of the numerator yields 



X [7 X 7].7'' + 27.[7 x (p - k)],7'^} . (21) 



Basic integrals over the loop momentum k can be easily evaluated to yield 
d'^k {l,k^} _ 1 



ivr^ A;2[(p — kY — rrfiY (1 ~ P)^ + P 



(22) 



where p = {m? — p^)/m^. Substituting these integrals into Eq. (j2T| and taking into account 
that [7x7] = 2z7o757 and {75,7^} = 0, we obtain 

S(p) = ;^^(P) (^^757.701^ - [7 X p].^ , (23) 



A(p) = -^fl + -^lnp). (24) 



where 

Ain\ = 

1-pV 1-p 

For s states, the angular integration in Eq. ()17j) is carried out by employing the following 
results for basic angular integrals (p = 1/2) 

X for /tj — X 

dpxU{p)^zXK,^{p) = { ' ' (25) 

1/3 , for K = 1 , 



— 2/3 i , for K = — 1 
2/3 i , for K = 1 , 



(26) 



where XKAt(p) is the spin-angular spinor 22| and cr denotes a vector consisting of the Pauh 
matrices. Finally, the result for the first vertex contribution reads (a is assumed to be an s 
state) 

a 



dprP^A{p) 



X <i ga{£aga + Prfa) - -fai^afa + PrQa] 



(27) 



where Pr = |p|, and = Qaipr) and /a = fa{Pr) are the upper and the lower components of 
the wave function, respectively. 

Now we turn to the second vertex contribution given by Eq. (fTBj) . The vertex function 
with two equal arguments Tji{p,p) can be obtained by the Ward identity. 



r(p,p) = -VpSW(p) 



(28) 



where T.^^^{p) is the free self-energy function introduced in Appendix El Simple differentia- 
tion yields 
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hip) 1 + b2{p) + h{p) p 
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(29) 
(30) 
(31) 
(32) 



In order to perform the integration over the angular variables in Eq. fjl8|) . we use the 
representation for the gradient in the spherical coordinates j3] 



d 



'p = PT^ + -Vn- 

OPr Pr 



(33) 



Now the angular integration can be carried out by using the following results for basic 
angular integrals (p = 1/2) 



dpxiuip) [P X ^n]zX^^^{p) 



, for K, = —1 
2/3 z , for fi; = 1 , 



(34) 



-4/3 i, for K = — 1 , 
, for K = 1 . 

The final result for the second vertex contribution reads (a is an s state) 

dpr pI 



(35) 
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X h{p){- 9afa + Qa'j^ - fa^ 

apr dpr 



-h{p){eafa + Pr9a)fa + h{p)fafa 



(36) 



This concludes our consideration of the zero-potential vertex contribution. 

The described approach to the evaluation of the zero-potential vertex term was first 
employed in P| . We mention also a different treatment of this term suggested by the Swedish 
group 0, 0|- In that work, the 5 function in Eq. (fT^ was approximated by a continuous 
Gaussian function with a small cutoff parameter. An advantage of the presented approach is 
that we end up with a single integration that remains to be performed numerically, while the 
consideration of P, 0] leaves a 4-dimensional integration in final formulas. This complication 
is not crucial for the zero-potential term since it is relatively simple. For the one-potential 
term, however, the consideration is more difficult, and the approach based on the integration 
by parts turns out to be extremely helpful. 

The zero-potential reducible term is given by 



(0) 
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where Qd is the lowest-order value of the (7-f actor given by Eq. 0. The derivative of the 
free self-energy function with respect to the energy argument reads 



de 
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(39) 
(40) 
(41) 



Integration over the angular variables yields 



^£1 = 3d{-^^] I dprpl 



oo 



X I ai (p) Eaigl + fa) + 2prgafa 

+a2{p){9l + f!)+as{p){gl-f!)y (42) 
Finally, the total zero-potential term Agi^^ is given by the sum of Eqs. (|2Z|), and (j^ . 

B. One-potential term 

The one-potential vertex contribution to the (7-factor can be written as 

X [A{p,r,p') X Vr6\p - p' - r)]ja{p') , (43) 

where the angular-momentum projection of the valence state is assumed to be = 1/2, 
p = {Sa, p), r = {Sa, t), p' = {Sa, p')? the function A is defined by 



A{p,r,p') = —Ania J 



d^k ^ — It + m 



(44) 



(27r)4 "'{p-ky-m? 

'6 — k — f + m '6' — k + m 

X — ^ — 7 

{p — k — r)2 — [p' — kY — ra^ 



and Vc'(r) = — 47ra;Z/r^ is the Coulomb potential in momentum space. Obtaining Eq. 
we have employed the vector potential in a form equivalent to Eq. ()14|1 

A(p - p' - r) = -'-{2nf [B x V.5\p - p' - r)] . (45) 



We mention also a factor of 2 included into Eq. (j4Hjl accounting for two equivalent diagrams. 
Integrating by parts, we divide the expression into two parts, 

= A^/lei,! + ^9^2 , (46) 

where 

A^?12,i = 4^m j ^^c(q)V^Jp)e.(p,p')^a(p'), (47) 



x[A(p,?,p')xS(q)],Va(p'), 



where g = p - p' = (0, p - p'), S(q) = VqVc(q), and 

e(p,p') = [Vr X A(p,r,p')]r=q " 



(48) 



(49) 



1. Agi^2 1 contribution 



We start our consideration of the Ag^^^,! term with the function 0, 

^ , ^ . f d^k ^- I/: + m 
0(p,p) = Ama I jt;^^!^— — j-^ ^7o 
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(p' — A;)^ — ' 



(50) 



where the gradient acts on the expression in the square brackets only. By using the following 
identity 



^' — 1^ + m f)' — 1^ + m ^' — 1^ + m 

^ (p' — kY — m? (p' — kY — m?'^ {p' — kY — rn? 



and commutation relations for the 7 matrices, we obtain 

d^k 



@{p,p') = —2ma J 



where the numerator q is 

Q 



k"^ [(p — kY — Tn?\ [(p' — kY — rn?\ 

lc{i> ^)7o{ (;^' - ^ + m) [7 X 7] 

+ [7x7] (/ -]ji^m)^-i" . 
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(51) 
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Now we employ the identity [7x7] = 2i7o757 and use the commutation relations in 
order to bring the matrices 70, 75, 7 to the right. This yields 



Q = 4i7o-(;^ - ^ + m) 

"^7 + (Po - ^o)7o7 - (P' - k) 



X 



7 75 



(54) 



Carrying out the summation over the repeating indices in the numerator, we obtain 



Q = 8i75 m^7 + 2mq - (p' - k)(;^ - ^) 

+ {p'o - kohloi^ - ^) 



(55) 



The integration over the loop momentum k is carried out by using the following results 

for basic integrals 

d'^k {1, /c^, kf^ki^y 
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i 1, xh^, X b^K - -Y9ny 



(56) 



where h — yp -\- {1 — y)p', N = xh^ + ym^p + (1 — y)'w?p' , p' = (m^ — p'^)/m^. 

By an explicit calculation one can show that the part of the basic integrals proportional 
to g^,^ yields a vanishing contribution. We therefore have 



@{p,p') = — / dxdy 
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m^-f + 2mq - Q'Q + QoTTo^?] , 



(57) 



where Q — p — xb — {1 — xy)p — x{l — y)p', Q' — p' — xb — —xyp + (1 — x + xy)p' . We now 
use the commutation relations in order to bring ]() to the left of 7 and to the right of 7. 
This yields 
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75^4o7 + (ClJ^ + C2;^')P 



+ {Dii) + D2i)')v + Fi iriol + F2 770/ 
+G'ip + G'2p' + ii"iP7o}, 



(58) 



where the coefficient functions are given by 



+ 2el{l - x)(l - xy) 
xy{l - xy) , 
- y) , 

— (1 — X + xy){l — xy) , 
(1 — x + xy)x{l — y) , 

- x){l - xy) , 
-£^(1 - x)x{l - y) , 
2m , 
—2m , 

-26^(1 - - xy) . 



(59) 
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In order to perform the angular integration in Eq. (|47jl . we define the set of scalar functions 
ri{i = l...Q): 

dxdy jpi X+{v)f^X+{v') + ^2 xl(p)crX-(p ) 

+P ^3 (P)X- (P') + P ^4 xl (P)X+ (P') 

+p'nxUp)X-(p') + P'nxUp)x+(p')| , (60) 

where x±(p) = X±fta,ma(p)- The functions Vi depend on = |p|, v'r = Ip'I; and ^ = p ■ p' 
only. They are given by 
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Cl{ea9a + P./a)/: + C2ga{eafa + K^?!) + {H, - G,)gX 
GliSafa + Prga)g'a + C^2/a(£a'7i + K/a) + (^1 + Gi)fag'a 
Di{eaga+Prfa)fa + D2ga{£afa+PWa) " <^25a/a 
I5l(£a/a + Prga)ga + D2fa{eag'a + p'rf'a) + G2fag'a 



(61) 



where Qa = ga{Pr), fa = fa{Pr), 9a = 9a{Pr) , f'a = fa{Pr) ^re the radial components of the 
valence wave function. 

Now the angular integration in Eq. ()47|1 can be performed, taking into account the fol- 
lowing results for basic angular integrals (/i = 1/2) 

^ { 271 f d^F , ioT K = -1 

/ dpdp'FxU{p)(T,x.,{p') = { -^-y (62) 

y -27r/3 y diiF, for k = 1 
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1, 



[dpdp'FxUp)p'zX-.,{p') = { {-i' ' ' (64) 

y -27r/3 j d^^F, for k = 1, 

where F is a function that depends on pr, p'^, and ^ only. Finally, we obtain the expression 
for the first vertex contribution, 

.2v roo ri ri ^2^/2 



^gier,! = J dprdp'^ J d^ J dx dy 
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2. Ag'^gj. 2 contribution 



First, we evaluate the function A in Eq. p8|) . Its expression can be obtained from the 
vertex function by the following identity 

A{p,q,p') = Vp,r\p,p'), (66) 

where T^{p,p') is defined Appendix 13 Using Eqs. ()A8j) - ()Alljl and taking into account that 

Vp.iV = 2(1-2/) i-xy p + (1 - a; + xy)p'] , (67) 
Vp/a = -xVp/A^ + 2(l-x)p, (68) 



we obtain 



A(p,g,p') = rfx dy I ^'-^^^ (^-xy p + (1 - x + xy)p'^ 
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+ ^ (^-2(1 - x)p7o + £aC7 + Diria^^ | 



(69) 



with the coefficient functions a and B-D given by Eqs. ()A9|) - ()A10|) . Note that the expression 
in the square brackets is similar to the corresponding term in Eq. ()A8|) . We can therefore 
write 
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where the expressions for jFi 2 can be obtained from ()A13j) . ()A14j) by a straightforward 
substitution. Finally, we obtain 



^a(p) A(p,g,p')V'a(p') = j rfxd^/|7^lxUp)'^X-(p') +^2X-(p)o-x+(p') 

+(p7^3 + p'7e4)xUp) x+(p ) + (p7^5 + pX)xI(p) x-(p ) 



(71) 



where 
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= ^a + w? + 2xN + eaH)gag', 

+eaB{eaga + Vrfa)g'a + ^aCgaiSaQa + V'rfa) 

+D{6aga+Prfa){ea9a+Prfa), (73) 
+eaB{eafa + Pr^a)/a + SaC fai^Safa + v'rQ'a) 

+D{eafa+Pr9a){eafa+P'r9'a)- (74) 

We perform the angular integration in Eq. PH|) taking into account the identity 

V^Vc{q) = 8'KaZ^, (75) 

Qr 

and the following basic angular integrals (/x = 1/2) 
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where F is a function that depends on pr, p^, and ^ only. The final result for the second 
vertex contribution reads 
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3. Reducible part 



The one-potential reducible contribution is given by 

dp dp' 



(1) f dpdp 

A^red = 9dJ -J^VM 



(79) 



where V^{p,p') is the vertex function, p = (e^, p), p' = {Sa, p'), and is the Dirac value of 
the (yf-factor. An expression for this contribution can be obtained by a simple differentiation 
of formulas for r^(p,p') given in Appendix 1X1 



C. Many-potential term 



Expressions for the many-potential term can be obtained by the point-by-point subtrac- 
tion of the corresponding zero- and one-potential contributions from the unrenormalized 
expressions ^ and 0. In order to transform them to the form suitable for a numerical 
evaluation, we have to perform the summation over the magnetic substates and the integra- 
tion over the angular variables. This can be carried out by using standard techniques (see, 
e.g., j|29|). We present here only the results for the unrenormalized contributions. The final 
expressions can be obtained from them by the corresponding point-by-point subtractions. 

The vertex correction to the (^-factor is given by (a is an s state) 



oo 



1/2 1/2 1 

J2 jl J 



P{ni,n2) Rjiuj,an2nia) 



where u = 1 — iO, Rj is the generalized Slater integral (the corresponding expressions 
can be found in |2£ 



P(ni,n2) = n(/i/20)(-l)'^ 



, /(2ji + l)(2j2 + i: 




n(/i/2/3) = [1 + (— 1)'^+'^+'^]/2, (• ■ • ) and {■ ■ ■ } denote 3j and 6j symbols, respectively, and 

/■oo 

Rn,n, = / drr'[gnAr)fn,ir) + /„,(r)^7„,(r)] . (82) 
The reducible contribution reads (a is an s state) 

^« r , 

Alfred = 

^ (-l)>-J"+-^P(a, a) Rjjuj, anna) 



nJ 



We mention that both the vertex and the reducible contribution are infrared (IR) diver- 
gent. The divergence occurs when en-^ = £n2 = in the vertex term and En = £a in the 
reducible term. However, the sum of these terms can be shown to be convergent. In prac- 
tical calculations, we perform the u integration for the sum of the vertex and the reducible 



contributions. In that case, the integrand is a regular function for small values of uj and can 
be numerically integrated up to the desirable accuracy. 



IV. HIGH-PRECISION EVALUATION OF THE IRREDUCIBLE 

CONTRIBUTION 



In this section we describe adaptations of the approach developed in [32, |33|, |3J| . In this 
method, the renormalization scheme is completely performed in coordinate space, which 
provide a great simplification in the numerical evaluation. In that method the integration 
over the photon energy z is done last. The complex plane integration contour is divided in 
two parts js^. The contour around the bound state poles of the electron propagator gives 
the low-energy part, which reads, in Feynman gauge: 



ir,L 



-2a p 



7r 



dz / dx2 / dx\ '\\}\{x2)a^G{x2^x\^z^a^j^ba{x\ 



^ sin[(£a - ^)x2l] 

X2\ 



B4) 



where 5a is the first order wave function correction ©, is the energy and i\}a the wave- 
function of the bound state under consideration and P denotes a principal value integral. 
Here we use the Feynman gauge for the low-energy part instead of the Coulomb gauge in 
Ref. |24l] to be compatible with the gauge used in the vertex and reducible part calculations. 
Here the integrand is well enough behaved that the transformation of the low-energy part 

jj K 



in Coulomb gauge done in 35| to improve numerical accuracy is no longer required. The 



high energy parts is written as 



Ai?ir,H 



—a 
vri 



/ dz 


^ dX2 J 


/ch -i 





-bx2i 



X21 



—26m / dx tljl{x)(36a{x) 



(85) 



where h = —i [(e^ — zY + iSf^"^ , Re(6) > 0, and X21 = X2 — Xi. The index fi is summed 
from to 3. The contour Ch extends from —ioo to — and from + ze to +ioo, with 
the appropriate branch of b chosen in each case. Note that the sign change in and (jHK|l 
compared to Refs. 1^4, 0| comes from the opposite convention G{e) = {H — e)^^ used in 
those works. 

In this part we use the Pauli-Villars regularization and 



Sm{A) = - 



-ln(yl2) + - 
4 ^ ^ 8 



(86) 



Because we do not do a transformation to Coulomb gauge in ()84j) angular integration is now 
identical to the one done for the high-energy part (jH3j) as described in js^, and lead to the 
same angular coefficients. 

The high-energy part is split into two pieces A£'ir,HA and Ai^ir^HB- In AEit-^ua the Coulomb 
Green's function is expanded around the free one up to the one potential term, the wave- 
function is expanded around Xi ^ X2, which is the place from which singularities that lead 
to ultraviolet divergence arise. This part must thus contains the regularization term. The 
regularized integral is evaluated analytically. The high-energy remainder AEw hr, which is 
finite is evaluated numerically. The details of the method are described in Refs. 
Because the cancellations that occurs in the normal self-energy (the contribution is formally 
of order am, while the final result is of order am{Za)'^, leading to the loss of 9 significant 
figures at Z = 1), are not present here, very high accuracy can be reached even at low-Z. 
At Z = 1 one looses only 3.5 significant figures for Is state and 4 for 2s states. Convergence 
of the numerical evaluation of the integrals is checked by doing a sequence of numerical 
calculations with increasing number of integration points in all three integrations. In the 
case of Z = 1, we have observed that the value obtained with this method, provides a final 
answer in slight disagreement with the results from the Za expansion, even though the con- 
vergence of the numerical integration would lead to think that accuracy is larger than the 
observed difference. This loss of accuracy is probably due to the code used for the numerical 
evaluation of the Green function in the high-energy part Ai?ir,HB- Due to that fact, we do 
not use the method described for this section for Z = 1, and for Z = 2, we used the half-sum 
of the value obtained with the method of this section and from the one described in Sec. IIIII 
The different contributions to the irreducible part, evaluated by the method described in 
this section are presented in tables |l] and |n] for Is and 2s respectively. 



V. NUMERICAL EVALUATION AND RESULTS 



We start with reporting some details about our numerical evaluation. The calculation 
for Is and 2s states was performed for the point nucleus. In addition, for the ground state, 
we repeated our evaluation for the hollow-shell model of the nuclear- charge distribution and 
tabulated the difference as the nuclear-size effect Ag^s- 

The calculation of the irreducible part is quite straightforward. For the point nucleus, the 



TABLE I: Contributions to the Is irreducible correction in unit of 10~^ (ppm). Note the cancel- 
lation between the low energy part and the piece of the high-energy part containing the renormal- 
ization terms. Numbers in parenthesis represent error in the last digit (when not given, accuracy 
is better than 1 on the last digit) 
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15 
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-5391 
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23 
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TABLE II: Contributions to the 2s irreducible correction in unit of 10~^ (ppm). Note the cancel- 
lation between the low energy part and the piece of the high-energy part containing the renormal- 
ization terms. Numbers in parenthesis represent error in the last digit (when not given, accuracy 
is better than 1 on the last digit) 
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5419 
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2 
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3 
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4 


71857 
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5 
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-5364 
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4 
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47 
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16 
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41131523977 


47 


267572 
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5387 
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perturbed wave function \6a) can be found analytically in a closed form In calculations 
involving an extended nucleus, the perturbed wave function was evaluated numerically, 
stored on a grid and then obtained in an arbitrary point by interpolation. The calculation of 
a non-diagonal matrix element of the self-energy function was carried out by two independent 
numerical methods. The first one is based on an expansion of the bound-electron propagator 
in terms of interactions with the binding field l3l|l (for details of the numerical procedure see 
), whereas the second one is described in Sec. IIVI Numerical results obtained by both 
methods are in good agreement with each other. Point-nucleus results for the irreducible 
part are listed in Tables ITTTl and HVl Presented values are obtained the by the second method 
in all cases, except Z = 1 and 2 for the Is state. For Z = 1, we used the result obtained 
by the first method; for Z = 2, a half-sum of the values obtained by the two methods was 
employed. 

The evaluation of the zero-potential and one-potential contributions Agi^^ and Agil^ is 
relatively simple. The zero-potential contribution [given by Eqs. pTj) . p6|) . and (j42|l ] con- 
tains a single numerical integration that can easily be carried out to arbitrary precision. The 
one-potential contribution [Eqs. (jHS)), ()78|1 . and ()79j) ] contains a 5-dimensional integration. 
The integration over one of the Feynman parameters can be carried out analytically, which 
speeds up the calculation significantly. The evaluation of this contribution is rather similar 

n 

to that of the one-potential term to the first-order self-energy correction 29] . Point-nucleus 
results for the zero-potential and the one-potential contribution are listed correspondingly 
in the third and the forth column of Tables lllll and IIVI 

The evaluation of the many-potential contribution is the most difficult numerical part 
of the present investigation and mainly defines the total uncertainty of the results. Since 
the actual calculation of the many-potential term is carried out by taking the point-by- 
point difference of the unrenormalized, the free, and the one-potential contributions, one 
should be aware about large numerical cancellations that occur in this difference. Additional 
cancellations appear when the vertex term is added to the reducible contribution. The 
final formulas for the many-potential term contain an infinite summation over the angular- 
momentum parameter. (This expansion is often referred to as the partial-wave expansion.) 
In our approach, we chose the absolute value of the relativistic angular parameter k of 
intermediate electron states, |k| = j ' + 1/2, to be the expansion parameter. The summation 
was evaluated up to l^maxl = 25 — 35, and the tail of the expansion was estimated by a 



least-squares inverse-polynomial fitting. 

The general scheme of the evaluation of the many-potential contribution is similar to that 
of our previous investigation Q]. A new feature introduced in this work in the case of the 
2s state is that we do not introduce subtractions anymore that remove infrared divergences 
separately in the vertex and the reducible term. Instead, we evaluate the integration over 
uj for the sum of Eqs. ()80|) and ()83p. In this case, infrared singularities that are present 
in two parts of the integrand cancel each other and the total integrand is a smooth regular 
function. This modification of the numerical procedure allows us to avoid an additional loss 
of accuracy due to numerical cancellations. 

The following contour was used for the uj integration in the present work: [eq — ioo, Eq — 
iO] + \eq — iO, — iO] + [iO, Eq + iO] + \eq + iO, £0 + This contour is advantageous for the 
evaluation of self-energy corrections in the low-Z region since it avoids the appearance of pole 
terms that lead to additional numerical cancellations. The parameter Eq in the definition of 
the contour can be varied. In actual calculations its value was taken to be about ZaEa for 
low Z . 

The results of our numerical evaluation for the Is and 2s states are presented in Tables ITTTl 
and IIVI respectively. For the ground state, we list both the point-nucleus and the extended- 
nucleus result. It should be noted that the uncertainty specified for the latter result refers to 
the estimated numerical error only and does not include the nuclear-model dependence. We 
estimate the model dependence of the nuclear-size effect Af^Ns to be about 1%. Numerical 
values of the root-mean-square radii used in our evaluation coincide with those of tTil. For 

TIQ 

the ground state, we compare our values with the results of the previous evaluations p, 12I|- 
The results of ^ were obtained for the point nucleus, while the evaluation |3] was carried out 
employing the homogeneously- charged spherical model for the nuclear-charge distribution. 

Finally, we compare our numerical values with the analytical results based on the Za 
expansion and isolate the higher-order contribution F\,,a.{Za) that incorporates terms of 
order [Za]'^ and higher. 



A^SE = - [1 + {Zaf a2o + (Za)^Fh,o.(^«)] , (87) 



where the first term in the brackets is the known Schwinger correction, the second term 
a2o was derived first by Grotch j^^l for the Is state and later generalized to ns states by 



Shabaev et al. [18[, 



The higher-order function Fh.o.(^<^) for 1* and 2s states is plotted in Fig. El 



VI. CONCLUSION 



In this paper we have presented our evaluation of the one-loop self-energy correction to the 
electron (^-factor of Is and 2s states in H-like ions. As compared to the previous calculations 
of this correction for the Is state, an improvement of accuracy of about an order of magnitude 
has been achieved in the low-Z region. For the most interesting experimental cases, H-like 
carbon and oxygen, our calculation improved the accuracy of the theoretical prediction for 

rn 

the (yf-factor by a factor of 3 for carbon and by a factor of 2 for oxygen [17], which reduced 
the uncertainty of the e 
for the electron mass is jj, |l2| 



ectron-mass determination based on these values. The new value 

nn 



me = 0.000 548 579 909 29 (29) (8), (89) 

where the first uncertainty originates from the experimental value for the ratio of the elec- 
tronic Larmor precession frequency and the cyclotron frequency of the ion in the trap, and 
the second error comes from the theoretical value for the bound-electron ^f-factor. 
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APPENDIX A: FREE ONE-LOOP FUNCTIONS 



The free self-energy function in the Feynman gauge and in D dimensions is given by 



•2t 



d^k 1 ]l)-^ + m 



UV divergences in the above expression are regularized by working in D — A — 2e dimensions. 
The mass parameter n is introduced in order to keep the proper dimension of the interaction 
term in the Lagrangian. We separate the UV-finite part of the self-energy function S^'* as 
follows 



(Al) 



where the mass counterterm is given by 



aC, 3 - 2e 

dm — — m . 

47re 1 - 2e 



and 



C, = r(l + e)(47rr(£^ 
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In the limit e — > 0, the renormalized part of the self-energy function is given by 
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An 



2m {1 + In p 
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(A5) 



where p = {m? — p^)/m^. 

The free vertex function in the Feynman gauge and in D dimensions is written as 



X 



^ ~ It + m ^ ^' — ^ + m 



{p — ky 



m^ 



{p' — ky 



m^ 



The divergent part of the vertex function can be separated in the form 

rnp,p') = ^7'^ + r^(p,y)- 



(A6) 



(A7) 



Here we present an explicit expression only for the time component of the renormalized 
vertex function Fp omitting terms of order e and higher and assuming that p° = = £„, 



Ol I A. 

^RiP,p') = 7^ / dxdy — 



27r 



N 
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+eaB]() + SaCi) + D^^jl + Bali 



(A8) 



where 



A = a + m'^ - N \^- + xlnN j , 
B = 2{l-xy){l-x), 

C = 2{l-x + xy){l-x), (A9) 
D = -(1 -x), 
H = -4m(l - x) , 

a = xy{l — xy)p'^ + x{l — y){l — x + xy)p''^ 

-2{1 -xy){l- X + xy){p- p') , (AlO) 



N = X [yp + (1 — y)p']^ + ym^p + (1 — y)m'^p' ■ (All) 

The integration over one of the Feynman parameters in Eq. ()A8|1 can easily be carried out 
leading to an expression, equivalent to that in [2^. However, we prefer to keep the vertex 
function in a more compact form ()A8|) here. 

For carrying out angular integrations, we introduce the scalar functions J-'i^2 that depend 
on Pr = IpI, p'r = Ip'I, and ^ = p ■ p' only. 



a 



V^a(p)r?j(p,p')^a(p') = ^ dxdy^ 

^1 xi (p) x+ (p ) + ^2 X- (p) x-{p 



X 



(A12) 



J^l = {A + eaH)gag'^ + eaB{eaga+ Prfa)g'a 
+eaCga{eag'a+Prf'a) 
+D{eaga + Prfa){eag'a + p'rf'a) , 

^2 = {A-eaH)fX + eaB{eJa+Vrga)f'a 
+eaCfa{eaf'a+p'r9'a) 
+D{eafa + Pr9a){eafa + PWa) ) 



(A13) 



(A14) 



where x±(p) = X±«a,m<,(p), and Qa = ga{p), fa = fa{p), g'a = 9a{p'), and = are the 

radial components of the valence wave function. 
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TABLE III: The one-loop self-energy correction to the Is-electron ^r-factor for H-like ions. All 
values are absolute contributions to the g'-factor (1/a = 137.035 989 5) and presented in units of 
10~^ (ppm). Individual contributions are listed for the point nuclear model. A^ns denotes the 
nuclear-size correction calculated for the shell model of the nuclear-charge distribution. The labels 
"pnt." and "ext." refer to the point-nucleus and the extended-nucleus result, respectively. 

Z Agi,. A^i? Agi]) A^i?^'' A^se (pnt.) AgfNs Ag-SE (ext.) 



1 


1.52928 2320.77563 


0.50250 


0.03305 


2322, 


.84046(10) 


0.00000 


2322. 


.84046(10) 
















2322. 


.8404(9)" 


2 


5.20640 2316.00970 


1.55757 


0.13053 


2322, 


.90420(9) 


0.00000 


2322, 


.90420(9) 
















2322, 


,9040(9)" 


3 


10.52313 2309.28506 


2.91759 


0.28869 


2323, 


.01447(9) 


0.00000 


2323, 


.01447(9) 
















2323. 


.0140(9)" 


4 


17.21613 2300.99753 


4.45945 


0.50260 


2323, 


.17571(9) 


0.00000 


2323. 


.17571(9) 
















2323, 


.1751(9)" 


5 


25.10744 2291.41521 


6.10392 


0.76661 


2323, 


,39318(9) 


0.00000 


2323, 


,39318(9) 










2323, 


.42(5)'' 




2323, 


.3928(9)" 


6 


34.06467 2280.73799 


7.79535 


1.07460 


2323, 


.67261(9) 


0.00000 


2323. 


.67261(9) 
















2323, 


.6724(9)" 


8 


54.78171 2256.69788 


11.16571 


1.79701 


2324, 


.44230(9) 


-0.00001 


2324, 


.44229(9) 
















2324, 


,4421(10)" 


10 


78.74380 2229.82629 


14.34905 


2.61754 


2325, 


.53668(10) 


-0.00002 


2325, 


.53666(10) 










2325. 


.28^ 




2325, 


.5355(10)" 


12 


105.51169 2200.79830 


17.21623 


3.48376 


2327, 


.00998(12) 


-0.00005 


2327, 


.00993(12) 
















2327, 


,0103(12)" 


15 


150.23524 2154.28732 


20.77184 


4.75746 


2330, 


,05186(16) 


-0.00011 


2330, 


,05175(16) 










2329, 


.79'' 




2330, 


.051(1)" 


18 


199.76448 2105.29703 


23.34254 


5.85654 


2334, 


.26059(20) 


-0.00029 


2334, 


.26030(20) 
















2334. 


.262(2)" 


20 


235.17645 2071.71455 


24.49971 


6.41815 


2337, 


.80885(24) 


-0.00052 


2337, 


.80833(24) 



TABLE III: Is g-factor. (continued) 



Z Agiy. Ag'^^ Agil^ A^i?^-* AgsE (pnt.) A^ns A^se (ext.) 

2337.50* 2337.86(1)" 
24 311.21758 2003.24295 25.53187 6.95150 2346.9439(3) -0.0015 2346.9424(3) 

2346. 92(1)°^ 

30 437.19791 1899.42146 24.22809 5.76765 2366.6151(3) -0.0057 2366.6094(3) 

2366.77* 2366.59(1)" 
32 482.18213 1864.94188 23.14976 4.73024 2375.0040(4) -0.0088 2374.9952(4) 

2374.97(1)" 

40 676.59132 1729.56158 16.47453 -3.19027 2419.4372(5) -0.0349 2419.4023(5) 

2419.45* 2419.39(1)" 

50 952.87040 1569.36938 5.18971 -22.43986 2504.9896(7) -0.1615 2504.8281(7) 

2504.09* 2504.827(8)" 

54 1074.52868 1508.85111 0.51786 -33.12979 2550.768(2) -0.282 2550.486(2) 

2550.487(8)" 

60 1270.36538 1422.34049 -6.00306 -52.2051 2634.498(3) -0.610 2633.888(3) 

2634.54* 2633.895(9)" 

70 1638.25837 1290.32234 -14.07335 -90.9411 2823.566(5) -2.193 2821.373(5) 

2823.39* 2821.39(1)" 

80 2072.57316 1174.47096 -16.61904 -135.1054 3095.320(10) -6.913 3088.407(10) 

3095.34* 3088.46(2)" 

82 2169.55512 1153.32686 -16.29951 -144.1073 3162.475(12) -8.687 3153.788(12) 

3153.85(2)" 

90 2601.21377 1075.78727 -11.90297 -178.5734 3486.525(20) -22.32 3464.205(20) 

3486.56(3)" 3464.35(2)" 
3487.30* 

92 2722.17025 1058.21204 -9.99010 -186.3945 3583.998(20) -28.12 3555.878(20) 

3556.05(2)" 



" Ref. Q, Ref. 0. 



TABLE IV: Various contributions to the one-loop self-energy correction to the 2s-electron ^-factor 
for H-like ions for the point nuclear model. All values are absolute contributions to the ^'-factor 
(1/a = 137.035989 5) and presented in units of lO"^ (ppm). 







Z A^ir Ag^ Ag^ 


A^i^^ A^SE 


(pnt.) 


2 


1 


4905 2320 


7711 





4471 





1317 2322 


8404(3) 


4 


5 


0539 2315 


9889 


1 


3379 





5245 2322 


9051(4) 


6 


10 


1863 2309 


2305 


2 


4286 


1 


1729 2323 


0183(6) 


8 


16 


629 


2300 


886 


3 


600 


2 


070 


2323 


185(1) 


10 


24 


209 


2291 


216 


4 


778 


3 


210 


2323 


413(2) 


12 


32 


796 


2280 


417 


5 


909 


4 


585 


2323 


707(2) 


14 


42 


290 


2268 


639 


6 


956 


6 


188 


2324 


074(3) 


16 


52 


609 


2256 


005 


7 


892 


8 


014 


2324 


520(3) 


18 


63 


683 


2242 


618 


8 


696 


10 


056 


2325 


052(5) 


20 


75 


453 


2228 


563 


9 


351 


12 


307 


2325 


674(5) 


24 


100 


887 


2198 


735 


10 


177 


17 


427 


2327 


225(5) 


30 


143 


170 


2150 


542 


10 


109 


26 


584 


2330 


405(5) 


32 


158 


234 


2133 


739 


9 


728 


30 


024 


2331 


726(6) 


40 


222 


652 


2063 


747 


6 


429 


45 


708 


2338 


536(8) 


50 


311 


098 


1971 


870 


-1 


448 


69 


823 


2351 


343(9) 


54 


348 


579 


1934 


243 


-5 


641 


81 


003 


2358 


184(9) 


60 


406 


838 


1877 


219 


-12 


891 


99 


640 


2370 


807(9) 


70 


509 


241 


1781 


353 


-27 


042 


136 


596 


2400 


149(9) 


80 


619 


038 


1685 


423 


-42 


849 


183 


153 


2444 


765(9) 


82 


642 


093 


1666 


316 


-46 


101 


193 


936 


2456 


245(9) 


83 


653 


785 


1656 


779 


-47 


729 


199 


543 


2462 


378(9) 


90 


739 


302 


1590 


411 


-59 


021 


243 


372 


2514 


064(9) 


92 


765 


177 


1571 


607 


-62 


163 


257 


585 


2532 


207(9) 



=0= 



(a) 



(b) 



(c) 



FIG. 1: Feynman diagrams representing the self-energy correction to the bound-electron gi-factor. 

The double line indicates the bound electron propagator and the wave line that ends with a cross 
denotes the interaction with the external magnetic field. 
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FIG. 2: The potential expansion of the vertex diagram. The single line indicates the free-electron 

propagator and the dashed line denotes the interaction with the Coulomb field of the nucleus. The 
terms of the potential expansion are referred to as the zero-potential, one-potential, and many- 
potential contributions. 
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FIG. 3: The higher-order self-energy contribution Fh.oX-^^ct) for the Is and 2s electron y-factors 
in H-like ions. 



